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Abstract A new two dimensional non-perturbative code to compute accurate oscil- 
lation modes of rapidly rotating stars is presented. The 2D calculations fully take 
into account the centrifugal distorsion of the star while the non perturbative method 
includes the full influence of the Coriolis acceleration. This 2D non-perturbative 
code is used to study pulsational spectra of highly distorted evolved models of stars. 
2D models of stars are obtained by a self consistent method which distorts spher- 
ically averaged stellar models a posteriori. We are also able to compute gravito- 
acoustic modes for the first time in rapidly rotating stars. We present the dynamics 
of pulsation modes in such models, and show regularities in their frequency spectra. 



1 Introduction 

Fast rotation is known to have a strong impact on stellar structure and pulsations: 
centrifugal distortion breaks the thermal equilibrium in a star and induces transport 
of chemical elements and angular momentum through large scale currents and shear 
turbulence, the centrifugal force distorts the pulsations cavity and the Coriolis force 
modifies the dynamics of modes. For many stars such as, for instance, 8 Scuti stars 
and Be stars, rotation has to be taken into account in order to use seismology as a tool 
to probe their interior So far the impact of fast rotation on observational frequency 
spectra is poorly understood, and some questions still need to be answered. What 
are the modes corresponding to each observed frequency? Which cavity inside the 
star do they probe? How can the physics inside this cavity can be extracted? All 
these questions call for realistic modelling. 
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2 The model 
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2.1 The 2D structure 

At the present time there are no fully two-dimensional stellar evolution models. 
Some progress is being made in the context of the ESTER project (for Evolution 
STEllaire en Rotation H |), but such models are not yet available. An alternative is to 
calculate spherically averaged stellar evolution models and then to distort the models 
a posteriori using a self-consistent method. We present here such an approach based 
on a method established by |J2|. It consists of building the acoustic structure of a 
rotating star in hydrostatic equilibrium. Poisson's equation together with the curl 
of hydrostatic equilibrium equation are alternatively integrated, taking into account 
the centrifugal force. This is done through an iterative scheme which starts with 
spherical profiles for the density and, given the angular velocity, leads to 2D profiles 
for density and gravitational potential. The pressure and Fi profiles are then found 
by solving the hydrostatic equilibrium and the equation of state. 



2.2 2D non perturbative pulsations 

For the computation of pulsations we only require the hydrostatic and dynamical 
structure of the star and not the thermal structure. We compute the oscillation modes 
as the adiabatic response of the structure to small perturbations - i.e. of the density, 
pressure, gravitational potential and velocity field - using the Eulerian formalism. 

The code developed for this purpose is named ACOR -Adiabatic Computa- 
tions of Oscillations including Rotation-. It consists in solving the hydrodynam- 
ics equations perturbed by Eulerian fluctuations, by direct integration of the two- 
dimensional problem. The numerical method is based on a spectral multi-domain 
method which expands the angular dependence of pulsation modes in spherical har- 
monics series, and whose radial treatment is particularly well adapted to the be- 
haviour of equilibrium quantities in evolved models at the interface of convective 
and radiative regions, and at the stellar surface. The radial differentiation is made by 
means of a sophisticated finite difference method which is accurate up to the fifth 
order in terms of the radial resolution. 

This code has been validated by comparison with the results of 15] for polytropic 
models. The agreement between the two codes is found excellent. 



3 Results 

We apply this modelling process to a 2 M,7, non barytropic, differentially rotating 
star, with a surface angular velocity of 80%i2i., where 12^ is the Keplerian critical 
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Fig. 1 Frequency vs Energy diagram of the pulsation modes computed for a 2D model of 2 Mq 
evolved star. 



angular velocity (12^ = -^/GM/R^, is the equatorial radius). This model star has 
been evolved from an initial composition of = 0.72 to = 0.35. The angular 
velocity varies with radius £2{r) with X2center = 3i2surface- Depending on the region 
of the frequency spectrum explored by our computations, five classes are found. 
The frequency vs energy distribution of these five classes of modes are presented in 
FiglH The kinetic energy is the Eulerian equivalent of the inertia of modes, and is 
given by: Ec{mode) = ff,dm \ v'{r) p / | v'(r = R) p. As in the non rotating case, 
modes trapped in the core have higher inertia than those trapped in the envelope. 
Therefore, this diagram gives access to the nature - i.e. gravity, pressure or gravito- 
acoustic- of the computed modes. 

At low frequency, represented by grey dots in FigJT] we mostly find modes 
trapped in the region of sharp variation of the mean molecular weight above the 
convective core, characterised by a sharp feature in the Brunt-Vaisala frequency. 
These modes are counterparts of high order g-modes in the rapidly rotating case. 
They have negligible amplitude in outer layers, and we expect them to undergo ra- 
diative damping, therefore, they should not be detected observationally. However 
the spectrum is overcrowded by this type of mode at low frequency, which makes 
the detection of modes trapped in the envelope in this frequency range very difficult 
in synthetic spectra. 

At higher frequency -green in Fig{T]-, we find whispering gallery modes, chaotic 
modes and island modes such as those computed in homogeneous models of stars 
by |i3J. Whispering gallery modes are the counterpart of high degree acoustic modes 
in the rapidly rotating case. They probe the outer layers of the star, but as they show 
very low visibility factors, they might not be detected. 

Chaotic modes -also green in FigH]- have no counterparts in the non rotating 
case. They are found only in very fast rotating models of stars. Their spatial dis- 
tribution does not present a simple symmetry. They have significant amplitude in 
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Fig. 2 Left: Echelle diagram for axisymmetric island modes. Right: Kinetic energy distribution 
of an island mode in a meridional plane. 



the whole stellar interior. The lack of symmetry should result in a low cancellation 
factor Therefore, these modes are expected to be detected observationally. 

Island modes -black dots in FiglT]- are the counterparts of low degree acous- 
tic modes in the rapidly rotating case. They probe the outer layers of the star, and 
present good visibility factors. Therefore, they should be easily detected observa- 
tionally. They are characterised by two quantum numbers i and h -the number of 
nodes along and perpendicular to the ray path, see Fig. |2]-. 

Some island modes -the most energetic island modes in FiglT]right panel- are 
gravito-acoustic modes -mixed modes- such as for example the £ = and the l~\ 
island modes which frequency is around 6.5X2^.. These modes have significant am- 
plitude in the stellar core as well as in the envelope. Therefore these modes hold 
very interesting seismic diagnostics, they should have a good visibility factor and 
carry information about both the core and the envelope. These modes could pro- 
vide valuable knowledge on the internal physics of rapidly rotating evolved stars 
including their differential rotation profile. 

In Fig. |2] (left) is given the Echelle diagram of axisymmetric modes computed 
for the model mentionned earlier This echelle diagram is very similar to echelle 
diagrams of acoustic modes in non-rotating stars, but it is built on different quantum 
numbers (see right panel of Fig. |2]i. This illustrates that with the appropriate defini- 
tons for the quantum numbers island modes spectrum is regular. This new kind of 
large separation should then provide a seismic diagnostic which is valid for rapid 
rotators H. Moreover, in the low frequency regime, the spectrum is less regular, 
and the ridges adopt a particular curvature. In the non rotating case, this curvature 
is attributed to the mixture of acoustic and gravity modes. In this case, where the 
star is at the middle of its life on the main sequence and is fastly rotating, one can 
wonder how rotation impacts this deviation from regularity. 
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